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ABSTRACT 

The majority of Milky Way extrasolar planets likely reside within a few kpc of the Galactic 
centre. The Galactic tidal forces acting on planets scale inversely with radius in the Galaxy 
and so are much greater in the inner Galaxy than in the Solar neighbourhood. Within a range 
of 3.5 to 10 kpc, the vertical tide from the Galactic disc is predominant. Interior to 3.5 kpc, 
the effects of the Galactic bulge cannot be neglected and the in-plane tidal components are 
as important as the vertical ones. Here, we quantify the orbital changes induced by these 
tides. We find that the greatest perturbations occur when the planetary orbit is severely 
misaligned to the parent star's orbit. When both planes are perpendicular, the eccentricity 
of the planet is driven to unity, although the semimajor axis is secularly unaffected. When 
both planes are coincident, the effect from Galactic tides is minimized, but remains non-zero. 
In these cases, we provide estimates for the survival times, as well as the minimum baseline 
eccentricity variation for all Milky Way exoplanets as a function of Galactic parameters. 
Inclinations similar to the Solar System's (« 60°) can easily cause eccentric Neptunes (at 
w 30 AU) around host stars deep within the Galactic bulge (within 50 pc) to experience 
eccentricity variations of several tenths, and cause the exoplanets with the widest-known 
separations (« 10 3 AU) to experience similar variations in the Galactic disc. These variations 
occur on timescales of a few Gyr, a fraction of a typical main sequence lifetime. 

Key words: planets and satellites: dynamical evolution and stability - planet-star interac- 
tions - fundamental parameters; The Galaxy: kinematics and dynamics - structure - disc 



1 INTRODUCTION 

The vast majority of the thousands of candidate and confirmed 
exoplanets reside in the Solar neighborhood Q0, which is ~ 8 
kpc from the Galactic centre. Given their abundance locally, 
it is natural to conclude that the whole Galaxy is teeming 
with exoplanets, in accord with the Copernican Principle. Some 
confirmation is provided by microlensing surveys, which typ- 
ically monitor source stars in the Galactic bulge. Intervening 
host stars can act as gravitational lenses, whilst their asso- 
ciated exoplanets can be det ectable as perturbations of the 
microlensing light curves (e.g. Dong et all 120091 ; Ijanczak et alj 
l20ld : lMivake et al.ll20"Tll : lYee et al.ll2012l ). In this way, exoplan- 
ets have been discovered at Galactocentric radii from ~ 3 to 
6 kp c, as exemplified by OGLE 2 007-BLG-050 dBatista et al.l 
120091 ) and OGLE-2003-BLG-235 (|Bennett et aJl hood ). Addi- 
tionally, the SWEEPS (Sagittarius Window Eclipsing Extra- 
solar Planet Search) collaboration identified 15 transiting ex- 
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oplanet candidates in the Galactic bulge, and concluded plan- 
ets ar e just as comm on there as in the Solar neighbour- 
hood (|Sahuet al .1120061 ). 

The density of the Galactic disc increases moving towards 
the centre in a roughly expone ntial manner with a scale le ngth of 
between 2 and 3 kpc (see e.g.. lBinnev fc Merrifieldll 19981 ). Addi- 
tionally, the innermost parts of the Milky Way are dominated by 
a 1O 1O M0 boxy-shaped bulge, as seen most promi nently in the 
COB E/DIRBE near-infrared light distributions (|Dwek et al.l 
|l995j). As the star density increases substantially towards the 
centre, there are many more possible exoplanet hosts in the inner 
Galaxy than in the remoter outskirts, like the Solar neighbour- 
hood. However, the Galactic environment itself becomes more 
extreme, as stellar collisions, encounters, and flybys are more 
frequent in the inner parts and the Galactic tides are stronger. 

The largest component of the Galactic tidal field in the 
Galactic disc acts perpendicular to the Galactic plane. Unfortu- 
nately, observations have not yet been able to identify a typical 
planetary orbit inclination with respect to the Galactic plane 
due to a strong bias: the majority of all exopl anet ary candi- 
dates ha ve been discovered by the Kepler mission (jBorucki et al.l 
l2011al lbl; iBatalha et al.l 12012). which observes a fixed patch of 
sky along the Orion arm and can detect candidates with only 
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nearly edge-on orbits. Programs to detect transiting planets 
see them edge-on at a wide variety of declinations and right 
ascensions, depending on the instruments' line of sights (see 
IVeras fc M occkcl 20JJ for a more detailed discussion). The ra- 
dial velocity technique, which is responsible for the discovery 
of the majority of confirmed exoplanets, provides no informa- 
tion about the planetary orbital inclination with respect to the 
Galactic centre. Further, the Solar System's invariable plane and 
ecliptic are misaligne d with the Galactic plane at an angle of ap - 
proximately 60° fe.g. lHuang fc Wadell966l ; lDuncan et al.lll9871 ). 
Although the Galactic inclination distribution of exoplanets is 
unconstrained, high inclinations certainly exist, and may even 
be typical. 

Regardless of their orientation and location in the Galaxy, 
all Milky Way exoplanetary systems experience the Galactic 
tide. Within its Hill or Roche surface, a star's gravity domi- 
nates and an ex oplanet may survive the ef fects of Galactic tides 
unscathed (e.g.. iHeisler fc Tremainelll986l ). Exterior to the Hill 
surface, the Galactic tides are always important. For a IMq 
star in the Solar neighbourhood, the Hill surface has an extent 
~ 10 5 AU, which is a measure of the size of any exoplanetary 
system. At a Galactocentric radius of 500 pc, the Hill surface 
is an order of magnitude smaller with a typical extent of ~ 10 4 
AU. For the planets in the Solar system, the Galactic tide is 
not generally important. Although effects of tides on planets in 
the Solar neighbourhood has been considered before, studies of 
tides on exopla netary systems in th e bulge has been restricted 
to Oort clouds |Brasser et al.ll2010l ). 

A rough rule-of-thumb is that the precession timescale due 
to tides Ptidc ~ P c 2 xt/-Ppi! where P ex t is the orbital period of the 
host star in the Galaxy and P p i is the orbital period of the planet 
around the star. For Jupiter, this gives Ptide ~ 10 14 yr, well in 
excess of a Hubble t ime. However, wide-separation planets are 
now known (see e.g., Goldman et al.ll2010] ; lLuhman et alJliblll ; 

2500 AU. 
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lj) with semimajor axes up to 
For such wide-separation planets in the inner Galaxy, Ptidc ~ 
10 yr or less. In other words, there can be significant effects 
from Galactic tides over the age of such exoplanetary systems. 



1.1 Wide Orbit Planet Motivation 

Given the observationally-motivated emphasis in this work on 
wide-orbit planets , here we provide an updated summary of 
I Veras et alf l|2009l) , which describes the prospects for the for- 
mation, survival and detection of these bodies. 



1.1.1 Formation of Wide-Orbit Planets 

Here we briefly review mechanisms for generating planets at 
different distances from t heir host stars. The core accretion for- 
mation mechanism (e.g. I Pollack et al.l 1 19961 ) can readily form 
planets at several AU but has difficulty forming planets at tens 
of AU. Even Uranus and Neptune, at ~ 19 AU and ~ 30 
AU, require particularly favorable cir cumstances to have been 
formed in situ f rom co re accretion dLevison fc Stewartl lioOl] ; 



iThommes et alj l2002h . iDodson-Robinson et all l|2009h claim 



that 35 AU is a rough limit beyond which massive gas giant 
planets must form by an alternate mechanism, such as grav i- 
tation a l inst ability in the disc (e.g. ICameronlll978l ; lBosslll997f ); 
iBolevi (|2009l ) claims that this limit is approximate ly 100 AU . 
Disc instability easily forms planets at tens of AU (Boss 2003, 



1201 ll ). and m ay or may not be able t o form planets at a couple 
hundred AU (|Bossll2006l ; lBolevll2009l ). 

Beyond a few hundred AU, planets are highly un- 
likely to have been formed in situ. Instead, they were 
likely formed and subsequently scattered outward within the 
same system due to multi-planet gravitational instability (e.g. 
Rasio fc Fordlll996l;IWeidenschilling fc Marzarilll99fj;lLin fc Idal 



19971 ). IScharf fc Menoul (120091 ) and I Veras et al.1 (|2009l ) showed 



that planet-planet scattering amongst massive gas giants which 
could have been formed by core accretion (a < 35 AU) may 
generate a population of planets from 10 2 AU - 10 5 AU. How- 
ever, those planets with the widest orbits are unlikely to remain 
bound for several tens of Myr if they inter act with any sur- 
viving tight-orbit planets. iBolev et al.l (|2012l ) performed multi- 
planet scattering studies assuming one planet was formed by 
gravitational instability at a = 100 AU. This initial wide-orbit 
planet induced significant radial mixing, ultimately producing 
a population of planets with separations between several hun- 
dred and several thousand AU (see their Fig. 15). Alterna- 
tively, wide-orbit planets could represent captured free-floaters 
(|Perets fc Kouwenhovenll2012l ; IVarvoglis et al.ll2012l ) , a possibil- 
ity that will become increasingly plausible with additional stud- 
ies t hat support the pu rported vast free-floating planet popula- 
tion jSumi et al.ll201ll ). 



1.1.2 Survival of Wide-Orbit Planets 

Planets scattered to stable wide orbits could remain in those 
orbits if they remain unperturbed. One potential destabilizing 
perturbative source is a remnant inner planet which survives the 
scattering phase. If the wide-orbit planet is sufficiently eccentric, 
then repeated perturbations by the inner planet will cause an 
additional scattering event and instability. This add itional scat- 
terin g event may not take place for tens of Myr (|Veras et al.l 
120091 ). 

Another perturbative source is the birth cluster itself. Al- 
though stellar flybys can actually aid the passage of a plane t 
onto a wide orbit (e.g. iMalmberg et al.ll201ll ; lBolev et al.ll2012l ). 
alternatively flybys could also cause instability and ejection. 
lAdams et al.1 |2006l ) estimate that over 10 Myr, for clusters with 
between 100-1000 members, the typical impact parameter for 
two stars passing each other is about 700-4000 AU. Any planets 
at or beyond these distances are likely to be severely disrupted. 
Whether this disruption triggers ejection or just a shift to a dif- 
ferent wi de orbit is dependent o n the geometry of the encounter 
(see, e.g. lVeras fc Moeckell2012l ). Planets with tighter orbits also 
may e xperience instabil i ties w ith a wide range of outcomes; re- 
cently |Parter^Quan3 |2012l) performed 10 Myr cluster simu- 
lations with single planets on circular orbits at 5 and 30 AU, 
and found escape, orbit disruption, or quiescent evolution were 
all possible outcomes for individual stellar systems. 

Overall then planetary interaction in clusters is highly 
model-dependent and disrupted planets may or may not sur- 
vive cluster evolution on wide orbits. A simpler argument for 
how planets can survive their birth cluster evolution is that in 
many cases the dissociation timescale for the cluster is much 
shorter than the timescale for a planet to experience a strong 
perturbation from a stellar fiyby, or even fully form. Although 
the cluster dissociation timescale is a func tion of both gas dy - 
namics as well as N-body dynamics (e.g. iMoeckel et al.l 20121 ). 



the dissociation timescale appears to be positively correlated 
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with the cluster membership population (e.g. ISterzik fc DuriserJ 
Il995l ; [Bate 20l3, and references therein). Therefore, the smaller 
the cluster, the more likely planets can survive the cluster phase 
undisturbed. 



We apply these equations in the regions dominated by the Galac- 
tic disc (Section^) and Galactic bulge (SectionJSJ) before charac- 
terising the minimum-possible eccentricity variation of exoplan- 
ets due to tides throughout the entire Milky Way in Section HJ] 
We discuss these results in Section [7] and conclude in Section [5] 



1.1.3 Detection of Wide-Orbit Planets 

Although the majority of exoplanets have been discovered with 
Doppler radial velocity variations or transit measurements, 
these techniques fail to probe the outer reaches of planetary sys- 
tems. The widest-orbit planets are instead discovered by direct 
high-contrast imaging. These detections often require follow-up 
imaging studies at different wavelengths and/or with different 
instruments to ensure that the star and companion are kine- 
matically associated with each other (have the same proper mo- 
tion) and to constrain the companion's mass. For example, the 
exoplanet with the widest-known orbit (a projected separation 
of ss 2500 Ain, WD 0806-661B b, was initial ly detected by 
iLuhman et al. 2011) but was later confirmed bv iLuhman et all 
l |2012ft . 



All 16 companions at separations of at least 100 AU from 
their parent star^f] have super- Jovian masses, and at least half of 
these are likely to be brown dwarfs by mass (|Spiegel et alj|201ll 
show the planet/brown-dwarf boundary to be 11 — 16Mj). For 
WD 0806-661B b, the companion's mass is reported as 8±2Mj, 
safely pinpointing this companion as a planet. Although the 
distinction between planetary masses and brown dwarf masses 
likely helps to indicate the way the objects are formed, the re- 
sulting difference in motion due to Galactic tides is negligible. 

Probing this super- Jovian mass regime is important 
for understanding the low-mass tail of the initial mass 
function and evaluating the extent of the brown dwarf 
desert (iMarcy fc Butler! l2000l; iGrether fc Lineweaverl 2006; 
iDeleuil et all l2008l ; iPinfield et al.l |2012| ). Further, exploring 
planetary systems beyond about 10 AU is important for compar- 
ison with the Solar System, and to asses s the prevalence of exo- 
Kuip er belts and exo-Oort clouds (e.g. iRavmond fc Armitagel 
l2012h . The Gemini Planet Imager (GPI), with first light due 
in January 2013, is dedicated to exploring outer regions of ex- 
osystems where radial velocity and transit surveys cannot reach. 
Forthcoming GPI discoveries will raise questions about the for- 
mation and fate of wide-orbit companions, and are likely to re- 
vise the current Galactic-wide exopla net population estimates 
l|Sumi et al.ll20lil ; ICassan et al.ll2012t h Our study helps to pre- 
empt this line of inquiry by considering the effect of Galactic 
tides on wide-orbit companions at all distances within the Solar 
Circle. 



1.2 Plan for Paper 

This paper examines how the Galactic tidal field may drive 
orbital evolution in exoplanets for disc and bulge host stars. 
We first rederive the equations of planetary motion subject to 
the Galactic tidal field in Section [2] We build a new three- 
component model of the Galaxy in Section [3] and demonstrate 
the relative importance of each component in different regions. 



3 Numbers are correct as of 18 November 2012 from the Extrasolar 
Planets Encyclopedia. 



2 THE GALACTIC TIDES 
2.1 Derivation 

The disturbances on a planetary orbit ca n be modeled with 
the pe rturbed two-body problem. Following iHeisler fc Tremainel 
(|1986T ). we use a non-rotating, rectangular coordinate system 
r = (x,y,z) centered on, and orbiting with, the host star. The 
star is assumed to move on a circular orbit of radius Ro about 
the Galactic centre with frequency £l G . Then the instantaneous 
unit vectors in the radial and tangential directions are 



= cos(fl G t) e x + sin(£l G t) e y 
= — sin(fiGt) e x + cos(fiGt) e.y 



(1) 
(2) 



The triad (e x ,e y ,e z ) has fixed spatial directions while 
(e x i , e y i , e z i ) rotates with the star. At the origin of both primed 
and unprimed coordinate systems, the gravitational force from 
the Galactic potential $ exactly balances the centripetal force 
for circular motion 



Q G Roe x , - V$ = 0. 



(3) 



At a general point in the exoplanetary system, the force per unit 
mass is the sum of the forces derived from the two-body problem 
of star and planet, the forces due to the rest of the Galaxy, and 
the fictitious forces caused by our choice of a non-inertial frame: 



F 



-e x , - V$ + n G Roe x > 



(4) 



Now, we Taylor expand the Galactic potential about the host 
star's position to obtain 



3>{R,z) 



$ [Ro 1 + 



2*' . x' 2 . y''^ 1/2 



R 



-ii-o 

12 

y 



tin tin 



* iRo+x ' + iR > z 



(3) 



This expansion implies further that 



<9$ ,J 
dR +X dR 2 J 



(Ro,0) 



9$ 



+ Ro {dR) (Ro 0) y 



+ 



<9$ | x ,(f$ 
dz dz 2 



e z , + 0(x ,y ,z ) (6) 

(fl ,0) 

If the host star lies in the Galactic plane, then the approximate 
symmetry z — > ~z of the Galactic potential ensures that the 
vertical gradient <9<E>/<9z vanishes. It is conventional to introduce 
the Oort constants 

'Rd£l G \ ^ 

(8) 



A = — 

' 2 dR 

b = - (n G + 

which gives us 
F 



(R ,0) 

Rdn G 
1 dR 



(«o.o) 



^e x , + (A-B)(3A + B)x'e x 



© XXXX RAS, MNRAS OOO.rTHMl 



4 Veras & Evans 



(A - B) 2 y'e y / 



(9) 



C«o,o) 



Finally, we can use Poisson's equation in cylindrical coordinates 
to relate the local density at the host star to the vertical gradi- 
ents: 



47rGp(i? ,0) -2(S 2 -A 2 



(10) 



g 2 $ \ 

9z2 J (iio.O) 

We th us obtain an equation derived by iHeisler fc Tremaind 



F = 



(A - B) 2 y'e y , - (4ttGp - 2(B 2 - A 2 )) z'e z 



(11) 



By converting from the primed to the unprimed coordinate sys- 
tem, this equation takes the form 



d 2 x 


G (m* + m p ) x 


~dP 


(x 2 +y 2 + z 2 f' 2 


d 2 y 


G (m* + m p ) y 


dt 2 


(x 2 +y 2 + z 2 f /2 


d 2 z 


G (m* + m p ) z 


dt 2 


{x 2 + y 2 + z 2 f /2 



~T~ K" 1 xx% ~h X xy%) 



+ T-y^X + T yyU 



+ r z 



(12) 
(13) 
(14) 



where m* and m p represent the masses of the star and planet, 
respectively, and the perturbations T are 



Txx = n 2 G {(l - S) cos {2fl G t) - 6} 

T xy = n%(l-6)sm(2n G t) 

T yx = fi G (l-5)sin(2fi G £) 

T yy = -nS[(l-<S)cos(2«ot) + <5] 

T zz = - UtvG P (R , 0) - 28Q G 2 ] 



(15) 
(16) 
(17) 
(18) 
(19) 



w here J = —(A — B)/(A + B). This form of the equations is given 
in crasser et al l | |201(M although our derivation makes it clear 
that the reference frame is not inertial, as they mistakenly claim. 
For any Galactic model, the Oort constants A and B, together 
with the shear 8, are calculable from the Galactic potential and 
vary with position in the Galaxy. 



2.2 Perturbative Equations of Motion 

The perturbed two-body problem is usually written down in 
terms of orbital elements, which provide greater intuition into 
how the size and shape of a Keplerian ellipse changes under the 
perturbations. We can express the equations of motion (|12[) - (|14[ 1 
in terms of orbital eleme nts by using Gauss' theory or Lagrange's 
equa tions (|Burnsl [l976l and pgs. 54-57 of iMurrav fc Dermottl 
1999). An alternative procedure for deriving the orbital equa- 
tions subje ct to very general pert urbative forces has been given 
recently by I Veras fc Evans! (|2012l ). 

A commonly-used approximation in non-linear systems is to 
separate the fast and slow oscillation variables. In celestial me- 
chanics, this procedure can be performed by averaging over the 
fast oscillations of the mean anomaly or true anomaly of the or- 
biting body. The result is averaged equations for the evolution 
of the more slowly oscillating variables, such as the semima- 
jor axis or eccentricity. This "adiabatic" approximation is valid 



There is a sign error in their expression for T z 



when Y/n 2 <C 1, where n is the mean motion of the planet. 
Effectively, for the Milky Way, this relation holds whenever the 
semimajor axis a < 10 4 AU. 

In the general case, a planetary orbit is not coplanar with 
the Galactic disc. In the Solar neighbourhood, the out-of-plane 
tidal component T zz is an order of magnitude greater than 
the in-plane components. Thus, the vertical tide usually dom- 
inates the motion, unless the inclination is very low. This ap- 
proximation has been used to simplify the equations in many 
previous studies IIHeisler fc Tr cmainc 1986; Matcsc & Whi tman! 
19891. ll992l;lMatese et al.lll995l ; lBreiter et al.lll996l ; lBrasser|[200ll ; 
Breiter fc Rataiczakll2005l ). The vertical adiabatic equations of 



= = 



motion are: 

' de s 
di i 

' di s 
It, 

<m 
lit 

du 
~~dT 



5eVT 



2n 

5e 2 sin 2uj sin 2i 



■ cos uj sin uj sin iT z 



8nVl - e 2 



cos i (2 + 3e 2 — 5e 2 cos 2lj) 



4nV 1 — e 2 
5 sin 2 u> (sin 2 i — e 2 ) 



-T, 



1-e* 



2nVl - e 2 



-T, 



(20) 
(21) 
(22) 
(23) 
(24) 



Planetary orbital elements appearing in these equations are the 
semimajor axis, a, eccentricity, e, inclination, i, longitude of 
ascending node, Q, and argument of pericentre, to. One may 
confirm th at these equat i ons reduce to other results found in th e 
literature (|Brasserll200ll : lFouchard||2004l : iFouchard et al.ll2006h . 
Equation (|20p is a consequence of the averaged Hamiltonian 
of the system being independent of the mean anomaly or true 
anomaly. 

Although these equations are useful for studies in the Solar 
neighbourhood, the planar tides are important for other loca- 
tions in the Galaxy. The planar adi abatic equations of m otion 
can be derived from the equations in IVeras fc Evans! l|2012l 'l and 
are written out in Appendix A. As tidal forces are derived from 
a gravitational potential, it is always true that T xy — T yx . This 
equality has the consequence that the planet's semimajor axis 
is never secularly affected by any tidal perturbations (see e.g., 
Eqs. l20l and lAl[l . Consequently, if the eccentricity tends towards 
unity, then concurrently the periastron tends towards zero and 
the apastron tends towards 2a. In this case, the planet will either 
collide with the star or leave the adiabatic regi me. The latter is 
likely to occur at distances of several 10 4 AU l|Veras fc Evans! 
2012), and must occur if the planet is to escape the system. 



3 THE MODEL OF THE GALAXY 

3.1 Galactic Bulge, Disc and Halo Models 

In all the numerical calculations in this paper, we use a three- 
component Galaxy model that reproduces observed local stellar 
kinematics data. The Galactic halo is represented by a logarith- 
mic potential of the form 

2 

$halo(iM = -f In (R 2 + z 2 q-' 2 + d 2 ) , (25) 

with Wo = 215 kms" 1 and d = 16 kpc (where R and z are 
cylindrical coordinates). The parameter q is the axis ratio of 
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the equipotentials, and controls whether the halo is spherical 
(q — 1), oblate (q < 1) or prolate (q > 1). The Galactic bulge is 
modelled as a Hernquist potential 

GM h 



Extent of The Hill Surfaces 



bulge 



(r) = 



r + e 



(26) 



using Mb = 3.6 x 10 10 M Q and e = 0.7 kpc. The bulge a nd halo 
are very similar to those used bv lFellhauer et all (2006). 

Many authors use Miyamoto-Nagai discs to represent the 
Galactic disc b ecause of the former' s simple form of th e po- 
tential (see e.g. iFellhauer et"ai1l2006l ; [Brasser et al.ll2010h . We 
have chosen not to do so here. Instead, we choose a more re- 
alistic exponential disc at the cost of more complex analytics; 
the quantities derived from the gravitational potential involve 
special functions. If the Galactic disc is razor-thin and exponen- 
tial with scale length Rd, its surface den sity and rotation curve 
are (Pg. 77 of lBinnev fc Tremaind 1 19871 1 



T,{R,z = 0) 



v cilc (R, z = 0) 



So exp 

ttGZqR 2 
Rd 



R_ 

' Rd 



io (A) A "° (2^ 



(27) 



(28) 



where Jo, Ii, Ko and Ki are modified Bessel functions. The 
scalelength Rd is 3 kpc, while the normalisation constant Eo 
is chosen s o as to reproduce the local column disc density of 
51M0pc -2 llFlvnn fe Fuchslll994l ). Although we have only given 
the res ult in the plane, the full three-dimensional potential is 
known (|Evans fc de Zeeuwlll992h . 

The superposition of these components gives a good repre- 
sentation of the Milky Way's rotation curve. The circular speed 
at the solar radius is ~ 220 kms -1 . The Oort constants at the 
solar radius are A — 14.5 kms -1 kpc -1 and B = —12.9 kms -1 
kpc -1 . These values are consistent with determinations of th e 
Oort constants from Hipparcos data l|Feast fc Whitelocklll997r ). 
which yield A = 14.8 ± 0.8 kms -1 kpc" 1 and B = -12A ± 0.6 
kms - kpc -1 . The three-component model therefore accurately 
reproduces all the local stellar kinematics within the quoted er- 
ror bars. Away from the Solar neighbourhood, there is consid- 
erable uncertainty in the values of the circular speed and the 
Oort constants. 

Lastly, the vertical tides depend on the density in the Galac- 
tic plane. The density of the bulge and halo can be straightfor- 
wardly generated through Poisson's equation. To generate the 
three-dimensional density of the disc, we take the column den- 
sity from Eq. (]28[) and smear uniformly ov er the scaleheight, h, 
of th e thin disc, where h w 300 pc (see e.g.. lBinnev fc Merrifieldl 
1998). Doing so is tantamount to assuming that the vertical disc 
distribution is exponential, for which there is good observational 
evidence. 

The Hill surface in the Galactic tidal field is estimated us- 
ing th e constancy of the Jacobi integral in lAntonov fc Latvshevl 
(|l972l ). To a good approximation, this is an ellipsoidal surface 
that is stretched along the line joining the star to the center of 
the Galaxy, and compressed in the two orthogonal directions, 
such that 



+ 



(29) 



The semiaxes in the Galactic plane are straightforward to 
find (jAntonov fc Latvshev|[l972T ) 




0.1 0.2 



8 12 20 35 



fl/kpc 



Figure 1. The semiaxes of the Hill surface in the x' (black), y' (blue) 
and z (orange) directions as a function of Galactocentric radius. No- 
tice that the Hill surface is most elongated in the x' direction (along 
the line to the Galactic centre) and most compressed in the z direction 
(perpendicular to the Galactic disc). 



( Cm, 



1/3 L , 2/Gm^ 1/a 

3 \ a 



(30) 



with a — 4A(B— A). The semiaxis in the direction perpendicular 
to the Galactic plane is more difficult to compute and an explicit 
formula does not seem to have been given before. We find: 



[Q(l+VT+Q)] 2/3 -Q 
[Q(l+yi+Q)] 1/3 



Gm, 



1/3 



(31) 



where Q = 4A(A — B)/T zz . 

The variation of the semiaxes of the Hill surface with Galac- 
tocentric radius are shown in Fig.[T] At the Solar neighbourhood, 
the surface has semiaxes (1.39, 0.93, 0.72) pc or (2.89, 1.92, 1.49) 
x 10 5 AU for a solar mass star. At 0.5 kpc from the Galac- 
tic centre, the semiaxes have shrunk by an order of magni- 
tude to (0.22,0.15,0.14) pc or (4.56,3.04,2.81) xlO 4 AU. The 
Hill surface is important as it gives the typical separation of a 
planet from its host star at which tidal effects become dominant. 
Within the Hill surface, a planet may nonetheless meander be- 
cause of tidal forces, but its wanderings will not usually lead to 
escape. Note too that the Hill surface is much more flattened in 
the Solar neighbourhood than in the Galactic bulge, because the 
matter distribution is more strongly dominated by the Galactic 
disc at the former location. 

Although our Galactic model is complicated, it is worth 
noting a helpful rule-of-thumb. Outside the inner kiloparsecs, 
the three-component model has an almost flat rotation curve 
with amplitude w 220 kms -1 . The circular frequency and the 
Oort constants are roughly given by 



A{R) w -B(R) 



110 kms" 



R 



Q G {R) = A{R) - B(R) 



220 kms" 
R~ 



(32) 



(33) 



Although we always use the full expressions derived from the 
potential of all three components in our numerical calculations, 
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Figure 2. A comparison of the amplitudes of the Galactic tidal forces 
throughout the Milky Way. The vertical tide is composed of three 
linearly additive components (bulge, planar and disc) all of which 
yield no net perturbation when i = 0°. The plot demonstrates that 
the vertical tide is at least 5 times as strong as the planar tide only 
in the range 3.5 kpc < R < 10 kpc. 



these simple expressions are useful in garnering physical intu- 
ition. 



3.2 The Galactic Regimes 

Here, we evaluate the contribution of each component in dif- 
ferent regions of the Milky Way. Doing so helps us to under- 
stand the dominant effects and to assess the validity of the 
popular practis e of neglecting the planar tides often used in 
previous works (Hcislcr & Trcmainc 1986; Matcsc & Whitman 
19891. 1 19921; fMatese et al.ll 19951 ; iBreiter et al.ll 19961 : lBrasserll200ll ; 
Breiter fc Rataiczakll2005h . 

We summarize the comparison in Fig. [2] We take pa = 
Pbuigo + Pdisc + Phaio and derive the three curves labeled "Bulge", 
"Disc" and "Halo" from this partition using Eq. (|19p with S — 0. 
The black "Planar" curve is traced from T = [Qg(R)] 2 , which 
represents the amplitude of the right-hand side in Eqs. (|15|l - 
(|18|l . The planar perturbations are a function of time, unlike 
the vertical ones. 

Figure [2] shows that the planar contribution to the total 
tide can reach at least 10% of the total in all regions of the 
Galaxy, and may reach « 50% where the bulge and halo are 
important. Previous studies' neglect of the planar tide is then 
justified at the ~ 10%-level, but only in the regime in which the 
tidal contribution of the Galactic disc dominates. This regime 
is typically over Galacotocentric radii between 3.5 and 10 kpc. 
Therefore, we use this approximation too in Section [4] for our 
study of the disc regime. Our bulge calculations in Section [S] 
however, all include the planar tide, in addition to the vertical 
tide. 

Depending on the accuracy sought, only one vertical com- 
ponent of the Milky Way (bulge or disc or halo) needs to be 
included for many regions of the Galaxy. However, around the 
bulge-disc transition region, at R ~ 1 kpc, and around the disc- 
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Figure 3. Planetary eccentricity evolution at different places in the 
Galactic disc for a Solar System-like inclination (i = 60°) and a wide- 
orbit planet (a = 1000 AU). In the upper panel, zuq = 0° for each 
curve. The lower panel demonstrates that the eccentricity may in- 
crease or decrease at different rates depending on the value of zdq. 



halo transition region, at about R 
must be included. 



20 kpc, two components 



4 THE GALACTIC DISC REGIME 
4.1 Fiducial Evolution 

The range 3.5 kpc < R < 10 kpc is the most straightforward one 
in which to explore exoplanet orbital evolution. In this regime, 
we need to account for only the contribution from the disc (pdisc) 
in the vertical tides (see Fig.[2]and Eqs. 1201124)) . The initial incli- 
nation of the planet can also be crucial in determining the effects 
of the perturbation. The vertical tide vanishes for i = 0° as T zz 
is vertically symmetric about the Galactic plane. However, there 
is little evidence to support coplanarity amongst planetary sys- 
tems and the Galactic disc. 

We then consider a planetary system with i = 60° - similar 
to the Solar System - with a one Solar-mass central star. We 
assume that the star evolves on the main sequence for about 10 
Gyr, a value which helps motivate the duration of our numeri- 
cal simulations. This duration is computationally achievable be- 
cause our simulations are adiabatic, and often feature wide-orbit 
single planets. 



© XXXX RAS, MNRAS OOO.riHTIl 



Exoplanets Beyond the Solar Neighbourhood 7 



The Polar Stationary Orbit Evolution 
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Figure 4. Characteristics of planetary systems highly inclined with respect to the Galactic plane, a possibly commonplace occurrence. For all 
plots, R = 4 kpc. For the top right plot and middle panels, ao = 2500 AU. For the bottom panels, ao = 5000 AU. The upper panel describes 
planetary evolution arising from polar orbits (io = i(t) = 90° or ig = i(t) = 270°) and a stationary argument of pericentre (ujq = ui{t) = oj cr ; t ). 
In this regime, given enough time, no planets survive. Analytical estimates for the survival timescale (Eq. I36H are shown as large colored dots at 
e = 1.0. The middle panels describe the motion when u)q is allowed to initially deviate from u> CT \ t . The result is still that no planets survive. The 
lower panels keeps wo fixed at w cr it, but allows io to deviate from 90°. The result is that planets can persist, although their eccentricities often 
approach unity. In the lowest panel, Qq = 0° is assumed. 



We provide representative eccentricity evolution profiles for 
wide-orbit planets residing at different locations in the Galactic 
disc in Fig. [3] The upper panel demonstrates the dependence 
on R with a fixed value of zuo; the lower panel illustrates that 
changing zuo can have a significant effect on the orbital evolu- 
tion. Over billions of years, bodies at a — 10 AU orbiting stars 
within the Solar Circle may become significantly more or less ec- 
centric. The smallest values of R typically cause greater pertur- 



bations, but not always; exoplanetary systems within the Solar 
Circle are likely to harbour more dynamically excited scattered 
discs and Oort clouds than the Solar System. Microlensing ob- 
servations of wide-orbit eccentric planets residing in the inner 
parts of the disc represent snapshots of dynamically changing 
systems. 
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Galactic Bulge Evolution 
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Figure 5. Planetary eccentricity evolution in the Galactic bulge. Each plot shows the evolution for a different (a, i) pair. The curves in both 
upper panel plots are coloured according to the legend in the upper right plot. All plots assume too = 0°, except the lower right plot, for which 
R = 50 pc and eo = 0.5 is assumed for all curves. This plot demonstrates that eccentric Neptunes deep within the bulge are dynamically active 
over the course of their main sequence lifetimes. 



4.2 Highly-inclined evolution 

As the inclination of a planetary system approaches 90°, the 
variation of the eccentricity is maximized. At i = 90°, the incli- 
nation does not chang e, and so polar orbits always remain polar 
orbits. Brass er] l|200ll ) found a stationary solution correspond- 
ing to ejerit = ± arcsin (s/TJE) ~ ±26.6°; this value may also be 
deduced from Eq. (|24p . The importance of this solution is that, 
whilst a, Q, i and cj do not change, the eccentricity increases 
and tends towards unity. Hence, this changing orbit potentially 
permits collision with the host star or may cause the planet to 
leave the adiabatic regime. 

Although a planet is unlikely to reside in this specific config- 
uration, it is useful to assess how nearby configurations evolve 
and how the planet's orbit is affected. We first show that a 
planet will eventually achieve e — > 1 for all ao and eo values 
in the top two panels of Fig. [4j given i(t) = io = 90° and 
oj(t) — ojo — arcsin (\/l/5). The resulting survival timescale is 
higher for tighter planetary orbits and more circular orbits. In 
the figure, we adopt R — 4 kpc, approximately halfway between 
the Earth and the Galactic centre. This choice is motivated by 
the microlensing planet searches, which are most sensitive to 
lenses half-way between observer and source star. 

Next, we consider small deviations from the stationary or- 



bit. We sample a wide variety of values of uuo, and show in 
the middle panels of Fig. [4] that for the initial values sam- 
pled, at i = 90°, luq asymptotically tends towards either o; cr it 
or (180° -l-cjcrit). Doing so eventually causes the eccentricity to 
tend to unity. Therefore, the value of ojo does not appear to 
affect the final outcome, just the survival timescale. 

By contrast, even if uiq — u cr it, deviations from i = 90° will 
prevent the eccentricity from reaching unity, but still periodi- 
cally increase its value to nearly unity. The bottom two panels 
of Fig. [4] illustrate this effect. In those panels, the planets sur- 
vive even if io = 85° . However, the eccentricities of those planets 
achieve values so close to unity that they may become unstable 
due to other factors, such as close passage to the parent star 
or a small impulsive kick from other Galactic phenomena. Note 
importantly from the bottom right panel that even the i = 50° 
planet, a full 40° astride from the stationary inclination solu- 
tion, exhibits eccentricity variations of several tenths. Further, 
the bottom left panel demonstrates that, although ui is not sta- 
tionary, it hovers around o; C rit, 180° — w C i-it or 180° + cj cr i t for 
the majority of the evolution. 

Returning to the strictly stationary orbit, we can estimate 
the survival timescale, t SUIV , analytically. By Taylor expanding 
the inclination about i = 90° to first order and the eccentricity 
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about e = to fourth order, we obtain an analytic solution to 
the first terms of Eqs. (|21|1 and (|22[) . The solution gives, for 
ujq = a; C rit and ujq — — cj C rit respectively, 



e(t) w ±e 



e(t) « ±e 



Veg + (2-eg)exp[2iT zz /n] 

\/2cxp [tT zz /n] 
v/2 - eg + eg cxp [2tT zz /n] 



(34) 



(35) 



In both cases, we obtain the same value of the survival timescale 
tsurv when setting e = 1, namely 



^sur 



2T 2 



In 



(36) 



Note that t SUIV « n/T zz at eo ~ 0.49. We overplot this ana- 
lytical estimate with large dots on the upper panel of Fig. [4] 
to demonstrate the quality of the approximation. As expected, 
the approximation worsens as eo tends towards unity instead of 
zero. 

Equation (|36[) suggests that any planet on an adiabatic po- 
lar stationary orbit has a finite survival time. For a main se- 
quence lifetime t m3 , the critical planetary orbital period around 
its parent star for which a planet does not survive is 



Tcrit 



^msT 



■hi 



ms 1 22 



2 - eg 



If we apply the disc density law from Eq. (|28[) . then 
h 



Tcrit 



111 



2 -eg 



2£msGIjO 

65400yr x In 
h 



e*p(R/R d ) 



2 -eg 



/ Cms \ 



0.3kpc 



cxp 



R - 8kpc 

Rd 



(37) 



(38) 



(39) 



For a star at R = 4 kpc, this yields a critical timescale on the 
order of 10 4 - 10 5 yr. 



5 THE GALACTIC BULGE REGIME 

Exoplanet evolution in the Galactic bulge regime is more com- 
plex due to the importance of both the planar and vertical tides. 
The vertical tide is composed of the contribution from the bulge, 
but for locations in the outer bulge (such as at R — 1 kpc), we 
must also include the contribution from the disc. Because of the 
interplay between the vertical and planar tides, attaining com- 
pact analytical results is difficult. Hence, we perform numerical 
simulations to explore phase space. 

The results are summarized in Fig. [5] All plots demonstrate 
a variety of eccentricity profiles; the upper plots model a wide- 
orbit planet (a = 10 3 AU) and the lower plots illustrate the 
evolution for a closer-in planet (at a = 100 AU and 30 AU). All 
plots assume vjq — 0°, except the lower right plot, for which 
R — 50 pc and eo = 0.5 is assumed for all curves. The difference 
in the evolution profiles from the upper plots (between the R = 
0.2 kpc and R — 0.3 kpc curves) arises from the initial sign of 
de/dt, which is determined by the initial relative magnitudes of 
Eqs. ([21]) and (A2). 

Wide-orbit planets are significantly affected at all locations 
in the bulge. In the innermost regions, the timescale of the eccen- 
tricity oscillations is fast: on the order of 1% of the typical main 



sequence lifetime. Also, the evolution in these innermost regions 
(within a couple hundred parsecs) becomes largely independent 
of inclination, as shown by the upper right plot. However, other 
locations in the bulge are affected by changes of inclination. The 
upper-right panel shows a similar evolutionary pathway to the 
upper-right panel of Fig. [4j even though here uj = 0° and not 

Wcrit- 

Tighter-orbit bodies, such as analogues of the trans- 
Neptunian object Ceto (a w 100 AU) and the planet Neptune 
itself (a ~ 30 AU), are affected significantly only if the host stars 
reside within the inner hundred parsecs of the bulge. The lower- 
left plot illustrates the drastic differences in eccentricity profiles 
due to shifting the value of R by just 150 pc. The lower-right 
plot presents the dependence on ujq for an eccentric Neptune at 
R = 50 pc. This planet cannot retain its primordial eccentricity. 
If, however, the planet was born on a more circular orbit, the 
extent of the eccentricity change would decrease. 



6 MINIMUM EXOPLANET ECCENTRICITY 

A planetary system whose invariable plane varies little from the 
Galactic plane may be modeled in the planar adiabatic limit. 
Assuming that the planet orbits the star in the same sense that 
the star orbits the Galactic centre ( "prograde" ) , the equations 
of motion become: 

da (40) 



dt 

de 5ef3;Wl - e 2 . r „ . n 

di 

Tit 



o 



(41) 
(42) 
(43) 



dm 5f2 G Vl - e 2 r , 

— = — ^— cos [2 (vj - Qct)} 

dt In 

We plot solutions of these equations in Fig. [6] for planets in 
the disc at R — 4 kpc (left panels) and in the bulge at R — 0.2 
kpc (right panels). We select eo = 0.5 for all plots except the 
lower-right plot. The upper panels show that planets in the disc 
regime typically feature eccentricity variations of 10 -6 to 10~ 5 
with a period of tens of Myr; planets in the bulge regime show 
variations of 10 -5 to 10 -4 with a period of a few Myr. These 
panels illustrate only eccentricity decreases because vjq — 0° 
for each of those curves. We sample wo at nine equally spaced 
values from and including 10° and 170° in the lower-left panel, 
demonstrating the sensitive dependence of the eccentricity evo- 
lution on vjq. This dependence is equivalent in the disc and 
bulge. The lower-right panel instead illustrates the dependence 
of amplitude on eo (for zuq = 0°). This relationship between 
amplitude and eo is not monotonic because of the eyl — e 2 
term in Eq. <|41[) . Consequently, the greatest variation occurs 
for e = l/y/2 w 0.7. 

This term prevents a complete analytical solution to Eqs. 
(|40[) - (|43| ). However, under the small eccentricity approximation, 
one may attain a closed analytical solution for the eccentricity 
evolution. The solution gives the following maximum eccentric- 
ity increase and decrease factors: 



1 ± 



5H G , x 9x KT 8 T . 2 
1 ± ^ kpc/yr 



(44) 



2n ' R 
where T = 2n/n is the planet's orbital period about its parent 
star and we used the approximations in Eqs. (|312[) - (|33[ ) . There- 
fore, all exoplanets in the galactic disc have eccentricities which 
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Figure 6. The minimum possible eccentricity evolution for a selection of exoplanets in the disc (R = 4 kpc) and bulge (R = 0.2 kpc) regimes. 
This limit is obtained when i = 0°; the evolution is according to Eqs. I !;:[ 1 ,'■'•[ . In the upper panels (with eo = 0.5), note the difference in both 
the frequency and amplitude of the oscillations. The initial semimajor axes for the curves in both upper panels are equivalent; in the lower panels, 
a = 1000 AU. The lower-left panel (also with eo = 0.5) demonstrates how the variations are phase-shifted according to -cuq. In the lower-right 
panel, otq = 0° is assumed for all curves, to demonstrate the dependence on eg; the amplitude peaks for eo = l/\/2 ~ 0.7. 



vary by a factor of at least (1 ± 5f2o/2n). Note that one cannot 
use Eq. (|44[) to determine when a planet will survive because 
this estimate is based on the adiabatic approximation and for 
low eccentricity. 

Now we consider how Galactic tides affect the pericentre 
advance or retreat of planets. The maximum variation of dvj/dt 
is equal to 



5Q 2 G 2 x 1(T 14 T 2 a 

— = - 2 radxkpc /yr 



(45) 



If we compute the critical semimajor axis at which the maximum 
pericentre precession rate from Galactic tides is comparable to 
that from general relativity, we obtain: 



(46) 



where c is the speed of light. We emphasize that these precession 
rate estimates are lower bounds because they were derived in the 
limiting case of planar adiabatic motion. 



7 DISCUSSION 

7.1 Implications from Tides 

Here, we discuss some of the implications of the tidal results 
that we have obtained from this work. First, the inclination 
of the planetary orbit with respect to the Galactic plane might 
be indicative of dynamical excitation in a planetary system, and 
vice- versa. Particularly, planetary orbits that are highly inclined 
to the Galactic plane will feature the greatest excitation. Al- 
though survival is likely for the smallest orbits in these systems, 
their eccentricity variations at a given (fixed) semimajor axis 
will be higher than in other systems. Despite this variation pe- 
riodically becoming zero, such periods of dynamical quiescence 
represent typically just a small fraction of the parent star's main 
sequence lifetime. Our lack of unbiased exoplanetary inclination 
data suggests that we cannot yet pinpoint a preferential plane- 
tary inclination with respect to the Galactic disc. However, our 
own Solar System and the variety of orientations exhibited by 
transiting planets prove that i ^ 40° can easily exist for plane- 
tary systems. 

Second, our results are not restricted to planets. Our anal- 
ysis may be extended to binary stars or belts of objects such as 
Kuiper belts or scattered discs. I Jiang fc Tremainel (|2010| ) study 
the evolution of wide binary stars in the Solar neighbourhood. 
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One can instead consider this evolution at other locations in the 
Galaxy by setting m p = m* in our Eqs. ()12p - ()14| ). For binary 
stars of equal masses, n will be increased by a factor of 2 1//2 , and 
hence, the time evolution of every variable (Eqs. 12111241 and [ Alt 
IA4[) will be decreased by a factor of 1.41. The amplitudes of the 
variations will otherwise remain unaffected. Modeling Kuiper 
belts and scattered discs involves imposing distributions of ini- 
tial conditions on our equations. Because the individual bodies 
in these belts are unlikely to interact with one another, they can 
be treated by our formalism. Oort clouds are typically too dis- 
tant to be treated in the adiabatic approximation. The ir evolu- 
tion must be modelled either with N-body simulations IjBrasserl 
l200ll ; iKaib et al.l 201lT) or w ith the nonadiabatic equations of 
motion (|Veras fc Evansll2012l ). 

Third, differential pericentre precession due to Galactic 
tides mi ght affect long-term N-body simulations of planetary 
systems. IVeras fc Fordl ((2010J) highlighted the danger of neglect- 
ing general relativity when modeling hierarchical multi-body 
systems with high relative inclinations, and the effect from 
Galactic tides might equally be important to incorporate, par- 
ticularly at high inclinations with respect to the Galactic disc. 

Fourth, improving observational precision might be able to 
place meaningful constraints on planetary systems. Currently, 
the smallest observational errors on planetary eccentricity mea- 
surem ents are about 1 x 10 -3 (|Wolszczanl Il994l ; IWelsh et al.l 
2012). Although these values are well above the baseline ec- 
centricity variations predicted by Eq. (|44l) . and are measured 
for tight orbits, they might be comparable to expected eccen- 
tricity variations for planet WD 0806-661, with a w 2500 AU 
l|Luhman et al.l l201lh. Two other planets with a > 1000 AU 
(e.g. iGoldman etaHboid : iKuzuhara et alj|201ll ) bolster theo- 
ries that planets can exist in such remote regions of planetary 
systems despite being formed elsewhere. 

Fifth, the influence of Galactic tides may represent another 
way in which planets can transform into Hot Jupiters. Near- 
polar planetary orbits might be highly eccentric (see Fig. [3). 
The pericentres of these planets might reside close enough to 
the star such that star-planet tides will damp both the semima- 
jor axis and eccentricity of the planet. The frequency of such 
transformed Hot Jupiters would be an increasing function of 
both age and Galactocentric distance. This mechanism is un- 
likely to be prevalent in the Solar Neighbourhood because the 
Galactic tidal timescales are likely to be too long to allow for 
star-planet tides to shrink the planetary orbits appreciably. 



7.2 Effects of Passing Stars 

Although stars collectively help establish the Galactic tide, in- 
dividually they can produce strong perturbations when they 
fly past a planetary system. The consequences of this brief 
perturbation may be comparable to tides acting over Gyr. 
The effect of flybys on planetar y systems has been studied 
before in the Galactic Disc (e.g. Zakamsk a fc Tremainel [20041 ; 
IVeras fc Moeckelll2012l ). Here we quantify the effect in the Galac- 
tic bulge. 

We use the formalism of IVeras fc Moeckell l|2012l ). which is 
sufficiently general to be applied to the bulge. They assumed 
every star has a single planet on an initially circular orbit, and 
computed cross sections for the planetary eccentricities to be 
perturbed by a given amount, k. We simplify their treatment 
by assuming each planet has a mass of lMj and initially resides 



Perturbations from Passing Stars 
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Figure 7. Eccentricity changes from passing stars in the bulge. Shown 
are the number of instances over a main-sequence lifetime that a 
planet's eccentricity should change by 10~~ 1 (orange lines) or 10~ 4 
(blue lines). Each set of four lines, from bottom to top (solid, dashed, 
dotted, dot-dashed) correspond to a samc = {1, 10, 100, 1000} AU. The 
plot shows that the eccentricity changes from stellar flybys in the 
bulge may be comparable to changes induced by the Galactic tides. 

at the same distance from their parent stars, a samc , all of which 
are have masses of IMq. Then the number of times, N, which 
a planet's eccentricity will be changed by k over a star's main 
sequence lifetime is equal to: 



AT = 



0905r? ( 



1000 AU 

(|Ae| > K,ij) 



Pbulg 



O.5M pc- 



10 10 yr 



(47) 



where the normalized cros s sections <T n orm can be read off di- 
rectly from Figs. 8-11 of IVeras fc Moeckell (|2012h . and n = 
Vdispersion/Ucriticai- We adopt a dispersion velocity of 200 km/s 
as representative of three-dimensional bulge dispersion veloci- 
ties. The critical velocity is the encounter velocity at which the 
total energy of the s ystem is zero and ionization is possible (e.g. 
iFregeau etai1l2004l ): 



= 2 



GMqMj 



(M + Mj) 



(48) 



We consider k = {l0" 4 ,10 _1 }, a same = {1,10,100,1000} 
AU and c ompute t m3 = 10.94 1 Gyr by assuming stellar metal- 
licity from lHurlev et all (J2000). Figure [7] demonstrates that the 
minimum possible eccentricity variations due to tides presented 
in Fig. [6] is likely to be surpassed by passing stars. In the oppo- 
site extreme, the eccentricity evolution of planets on wide polar 
orbits appears to have comparable contributions from tides and 
stellar flybys (see Fig. [5]). Therefore, both of these effects may 
be equally important and deserve detailed study. 



8 CONCLUSIONS 

The effects of Galactic tides have been studied before in the So- 
lar neighbourhood. This is the first investigation into the effect 
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of tides on planetary systems within the Oort cloud throughout 
the rest of the disc of our Galaxy, as well as in the bulge. 



In the Solar neighbourhood, the Hill surface has an extent of 
~ 10 5 (M/M o ) 1/3 AU. For the Sun, this marks the outer bound- 
ary of the Oort Cloud, and Galactic tides are known to play an 
impo rtant role in dislodging comets into the inner Solar Sys- 
tem ( Smoluchowski fe Torbettlll984l; iHeisler fe Tremaind 1 19861 ; 
Duncan et al. | 



19871; iMatese fe Whitman! 1 19891; iLevison et al.l 



20011 ; iKaib fc QuinnlT2009h . However, the effect of Galactic 
tides on the Sun's planets is negligible, even over the age of 
the Solar system. However, the architecture of exoplanets is 
very varied, and wide-orbit exoplanets with semim ajor axes be- 
tween 100-2500 AU are known to exist (see e.g. , iKalas et al 



20051 ; iGoldman et~ail20ld : iLuhman et a ]|201ll ; lKuzuhara et a 



201 If ) . In the inner Galaxy, the tidal forces are stronger, typi- 
cally scaling inversely with radius in the Galaxy. For example, in 
the inner 500 pc, the Hill surface has a size of ~ 10 4 (M/M Q ) 1/3 
AU, whilst the timescale over which evolutionary effects become 
noticeable is billions of years. 

The Galactic tide is dominated by the contribution from 
the Galactic disc over radii in the range 3.5 to 10 kpc. In this 
disc regime, the popular practice of neglecting the planar tides 
is justifiable. Planetary systems that are at least moderately 
inclined to the Galactic disc are susceptible to slow but signifi- 
cant planetary eccentricity evolution over the parent star's main 
sequence lifetime. At high inclinations to the Galactic plane, ec- 
centricity evolution is greater. As the semimajor axis is never 
secularly affected by tidal perturbations, the periastron tends to 
zero and the planet is driven towards the host star on timescales 
of > 10 9 yr. 

Within 3.5 kpc of the Galactic Centre, the contribution of 
the Galactic bulge is always important. As the matter distri- 
bution is more spherical, the planar tidal components are as 
important as the vertical ones. Evolution in the bulge regime 
is now a result of complex interplay between the vertical and 
planar tides. Wide orbit planets (a = 1000 AU) are substan- 
tially affected by tides on timescales (tens to hundreds of Myr) 
that are much smaller than the main sequence lifetime (~ 10 
Gyr). Eccentricity variations of several tenths are very typical 
for these planets. They are never in a state of quiescence; their 
orbital parameters are continually changing though the effects 
of tides. Closer-in planets, at e.g., 100 AU, are only affected if 
the host star resides within the inner hundred parsecs. 

Our study of tides is a first step in the understanding of how 
exoplanetary systems interact with their Galactic environment. 
The main limitation of our work is that the host star has been 
assumed to move in a circular orbit in the Galactic plane. In fact, 
stars lead much more exciting lives! They are usually inclined to 
the Galactic plane, their orbits are usually eccentric and some- 
times chaotic, and they suffer perturbations that can move them 
many kiloparsecs. There is, for example, both chemical and dy- 
namical evidence that the Sun has been moved substantially 
from its place of birth (|Clavtonlll997l ; ISellwood fc BinnevllaXM) 
by spiral waves. There are therefore good reasons for believing 
that the effects of tides are still more substantial than we have 
found here! 
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APPENDIX A: PLANAR ADIABATIC EQUATIONS 

We obtain th e plan ar adiabatic equations for a planet subjected to tidal perturbations T by setting T xy = T yx in Eqs. (25)-(29) 
of I Veras fc Evans! (|2012l) : 



'da- 
dt , 



(Al) 



v 



' de\ 5evT" 



[4 cos i cos 2uj sin 20 + sin 2ui cos 20 (3 + cos 2i)] (T xx + T yy ) + 2 sin 2 i(T xx — T y 



dt I 16n 

/ P 

— 2 [4 cos i cos 2co cos 20. — sin 2a; sin 2S7 (3 + cos 2i)] T xy ^ (A2) 
1 1 [sin 20 (2 + 3e 2 + 5e 2 cos 2a;)] (T„ — — 10e 2 cos £ sin 2a; sin 2 Sl(T n + Tj, 



e 



+ 20e 2 cos i sin 2a; cos f2 sin OT xy \ ( V-"!) 

[sin 2 f2 cos i (—2 — 3e 2 + 5e 2 cos 2a;)] (T xx + T yy ) + 5e 2 cos O sin f2 sin 2w(Y xx — T y 



'dO\ 1 
, dt J p ~ invT^P ^ 

— 2 [cos SI sin O cos i (— 2 — 3e 2 + 5e 2 cos 2a;)] T xy ^ (A4) 
' £7u(T M + T ra ) + (2C*io + Cn) (T M - T w ) - 2 (2C*io cot 20 + C 13 ) T xy \ (A5) 



\dt J p 16nVT^ 
where, using the notation of I Veras fc Evans] (|2012T ). we have 

do = 5 (e 2 - 2) sin (2a;) sin (2fi) cos i (A6) 

Cia = cos 2ft (l - 6e 2 - 5 (-3 + 2e 2 ) cos 2a; - 10 cos 2i sin 2 a;) (A7) 

C12 = ll-6e 2 +cos2a;(5-10e 2 ) + 10 cos 2i sin 2 a; (A8) 

Ci 3 =s sin2fi(-l + 6e 2 + 5(-3 + 2e 2 )cos2a; + 10cos2isin 2 a;) (A9) 
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